Comparing Criteria for Circular Orbits in General Relativity 
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We study a simple analytic solution to Einstein's field equations describing a thin spherical shell 
consisting of collisionless particles in circular orbit. We then apply two independent criteria for 
the identification of circular orbits, which have recently been used in the numerical construction of 
binary black hole solutions, and find that both yield equivalent results. Our calculation illustrates 
these two criteria in a particularly transparent framework and provides further evidence that the 
deviations found in those numerical binary black hole solutions are not caused by the different 
criteria for circular orbits. 
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Binary black holes are among the most promising 
sources of gravitational radiation for the new generation 
of gravitational wave detectors LIGO, VIRGO, GEO and 
TAMA. Motivated by the need of theoretical models for 
the identification and interpretation of future gravita- 
tional wave signals, several researchers have solved the 
constraint equations of Einstein's field equations to con- 
struct initial data describing binary black holes in quasi- 
circular orbit 0,|,|,|,|, | . 

Constructing such initial data requires making sev- 
eral choices, including the decomposition of the initial 
value problem and the background geometry and topol- 
ogy. Moreover, solving the constraint equations provides 
the gravitational fields for black holes with arbitrary sep- 
aration and momenta, and an additional criterion has 
to be applied to identify circular orbits. It is not sur- 
prising that different choices lead to physically different 
data. While all of these different data may be correct so- 
lutions to the constraint equations of general relativity, 
some may be more relevant astrophysically than others, 
in that they better represent a binary black hole system 
as it arises from inspiral from large separation. 

The results of Cook Q and Baumgarte (2| (which we 
will jointly refer to as CB) and Grandclement, Gourgoul- 
hon and Bonazolla [Q (hereafter GGB) differ by about a 
factor of two in the orbital frequency for the innermost 
stable circular orbit. This discrepancy raises two ques- 
tions, namely which results are more relevant astrophys- 
ically and which choice in the respective approaches are 
responsible for the deviations. The better agreement of 
the GGB results with post-Newtonian results suggests 
that these represent binary black holes in circular orbits 
more accurately There is also increasing evidence 
that the differences between CB and GGB are related to 
the different decompositions of the constraint equations 
JTT| , [T^ | . CB adopt the conformal transverse-traceless de- 
composition, which allows for an analytic solution of the 
momentum constraint jl3| , while GGB adopt the confor- 
mal thin-sandwich decomposition [l5| [l(| (see also 
p2|). It has been demonstrated that the two decompo- 
sitions may lead to physically different data (Tl), and it 
has also been suggested that the thin-sandwich decompo- 
sition together with maximal slicing may provide a more 



natural framework for constructing quasi-equilibrium so- 
lutions |l6| . 

In this Brief Report we explore the effect of another dif- 
ference in the approaches of CB and GGB, namely the 
criterion for locating circular orbits. CB adopt a turning- 
point method, in which circular orbits are identified with 
extrema of the binding energy (see eq. ( |27| ) below) , while 
GGB identify circular orbits by equating the ADM (rif| 
and Komar masses jl8| (eq. (|3l|)). Since the two mass 
definitions agree only for stationary spacetimes, this cri- 
terion is closely related to imposing a relativistic virial 
theorem [ fl9| . 

To explore the effect of these different criteria for circu- 
lar orbits in a particularly simple and transparent frame- 
work, we apply them to an analytic solution of Einstein's 
equations describing a thin, spherical shell of identical 
collisionless particles. At every point on the shell the 
particles move isotropically, but all with the same speed 
in the plane perpendicular to the radius. In an oscillating 
shell, each particle moves about the center in a bound or- 
bit. In the Newtonian limit, each orbit is a closed ellipse, 
and for static shells each orbit is circular (compare [M ) . 
Since each particle follows a geodesic, circular orbits can 
be identified without ambiguity. These orbits can then 
be compared with those obtained from the turning-point 
and mass methods. 

In the following we will focus on a moment of time 
symmetry, when at least momentarily each particle is in 
a purely tangential orbit u r = (where u a is the four- 
velocity). The spherically symmetric line element can 
then be written as 



ds 2 



-a 2 dt 2 + iP 4 (dr 2 + r 2 (d0 + sin 2 6d<j>)), (1) 



where a is the lapse function and ij) the conformal factor. 
The rest mass Mq of the shell can be computed from 



M 



p u 



4tt / p Wijj b r 2 dr, (2) 



where g is the determinant of the spacetime metric and 
where we have defined the particles' Lorentz factor W = 
—au*. Since the shell's co-moving density po, which is 
a sum of the individual particle densities Pq, vanishes 



2 



everywhere except at the radius R of the shell, we find 



»- Erf -5sw*» (3) 

The conformal factor tp in ([!]) can now be found from 
the Hamiltonian constraint 



M 



V 2 ip = -2ni{j 5 p N , 



(4) 



and, following GGB, the lapse a from the maximal slicing 
condition 

V 2 (aip) = 2naij 5 (p N + 2S). (5) 

Here pn is the density measured by a normal observer 
n a 

Pn = n a n b T ab = n a n b £ p$u£u£ = p Q W 2 , (6) 
A 

(compare & |2l[] ) , and S is the trace of the spatial stress 

S = 7^ = P oi ll u t u 3 = p Q (W 2 - 1), (7) 
where we have used the normalization condition 



1 



W 2 - -f lj UiUj. 



(8) 



For time symmetry both the momentum density j a = 
and the extrinsic curvature vanish, so that a 
zero shift ji % = identically satisfies the shift equation 
obtained in the conformal thin-sandwich decomposition. 

The Hamiltonian constraint @ and the maximal slic- 
ing condition (^|) can readily be solved analytically by 
matching two vacuum solutions at the shell's radius R. 
Choosing the vacuum solutions such that the interior so- 
lution is regular at the center, while the exterior solution 
is regular at infinity, we find for the conformal factor 



W 



W\rR 

w 



for < f < R 



for r > R. 



(9) 



Here and in the following we non-dimensionalize all quan- 
tities with respect to Mq, e.g. f = t/Mq. The value of 
ij)\ R can be found by evaluating the conformal factor at 
f = R, which yields a quadradic equation with the solu- 
tion 



1 I W 
^ = 2 + Vl + 2l- 



(10) 



The sign has been chosen so that ip approaches the grav- 
itational potential </>Newt in the Newtonian limit. 



In terms of the ADM mass 17 



M A 



1 



DM 



2nM 



D^d 2 S t = (11) 
V>\R 



the exterior conformal factor (||) can be written 
Madm 



2f 



for r > R. 



(12) 



The maximal slicing condition (^) can be solved anal- 
ogously to the Hamiltonian constraint, yielding 



aip 



1 -^W 2 -_2) iQiQ<f<R 
2W^\ R R ~ 



2Wtp\ s r 

Dividing by ip, we find in the exterior 

= -a\ R {iW 2 -2) + 2Wj ) \ R f 
a W 2 + 2W^\ R f 



for ¥>R. (14) 



Evaluating this expression at r = R determines the coef- 
ficient a\ R 



a s = 1 



2W - 1 



(15) 



Following GGB we now compute the Komar mass pi 



Mr 



47rAfn 



Wad 2 Si 



a\ R (SW 2 — 2) + W 2 



(16) 

We note that the Komar mass is a slicing dependent 
quantity, and that this particular form results from hav- 
ing imposed maximal slicing. In terms of the Komar and 
ADM masses, the exterior lapse a can be written 



2f - (2M K ~ Madm) 
2f + Madm 



for r > R. 



(17) 



This expression reduces to the lapse as identified from 
the Schwarzschild metric in isotropic coordinates (see, 
e.g., exercise 31.7 in ^3|) only if the two masses agree, 
Mk = Madm (compare criterion (|3l|) below). 

So far, the shell's radius R and Lorentz factor W ap- 
pear independently in the above equations. It is intu- 
itively clear that searching for circular orbits will yield a 
relation between the particles' angular velocity and the 
gravitational field, and hence between R and W. Since 
our model consists of collisionless particles, circular orbits 
can be determined directly by solving the geodesic equa- 
tions. Since all particles are identical, it is sufficient to 
evaluate the equation of motion for one particle, which 
we take to orbit in the equatorial plane. We therefore 
have u e — u r = 0, so that the normalization condition 
(H) yields a relation between and W 



(«*) s 



W 2 



^IrR 2 ' 



(18) 



where we temporarily drop the bar notation. We now 
evaluate the geodesic equation, 



du a 
~d\ 



T a bc u b u c = 0, 



(19) 



3 



for a = r to find a condition for the particles to remain 
in a purely tangential orbit (du r /dX = 0) 

r^) 2 + r W u*) 2 = o. (20) 

Combining this with (fL8ft and W — era* gives 



W 



a 2 T- 



(21) 



When evaluating the Christoffel symbols, we must take 
into account the discontinuity in the first derivative of 
the metric coefficients at r = R. By averaging such a 
quantity over an extended shell and letting the thickness 
of the shell go to zero, we find that the derivative has to 
be replaced with 



V\r -> 7jW , ,r| + + lp,r\-) = ^Vv| + - 

Using this rule for both ip and a we find 



M A 



DM 



DM 



and 



1 tt 



Ah 



1 



M A 



DM 



2R 



1 - 



R 



M K 



2R 2 V 2R J V -R + Madm/2 
at r = R. Inserting these into eq. (|2l]) yields 

M K 



W = 1 - 



2R - 2M K + M A 



DM 



(22) 



(23) 



(24) 



(25) 



After some algebraic manipulation and dividing out the 
unphysical root W = 0, eq. (25) can be expanded into 



AW 5 - 6RW 4 - AW 3 + 10RW 2 + W - AR = 0, (26) 

where we have reintroduced the bar notation. This is 
the condition relating W and R for circular orbits. It is 
easy to show that this equation reduces to f2 2 = R~ 3 /2 
in the Newtonian limit (with > 1, d < 1 and W ~ 
l+v 2 /2=l + R 2 ti 2 /2). 

For black holes, alternative criteria have to be used to 
identify circular orbits. In the following we will compare 
the turning-point method adopted by CB and the mass 
criterion adopted by GGB. 

In the turning-point method, a circular orbit is identi- 
fied by finding an extremum of the ADM mass (or equiv- 
alently the binding energy) at constant angular momen- 
tum U<j, 



dM A 



DM 



dR 



= 0. 



(27) 



In a Newtonian context, this condition arises natu- 
rally from Hamilton's equations of motion. We start 



by differentiating the normalization condition, (u^) 2 
ip 4 \ R M$R 2 (W 2 - 1), with respect to R to find 



dW _ -(W 2 -l){l + b) 
~d~R ~ RW(1 + b)+ AW 2 - 2 



(28) 



for sequences of constant angular momentum, where for 
convenience we have abbreviated b = (1 + 2W/R) 1 ' 2 . 
We now locate an extremum of the ADM mass ( |ll| ) by 
setting its derivative with respect to R equal to zero 



W 



dW 



dR \Rb(l + b) 



1 



W 



R 2 b{l + b)' 



Combining ([28|) and ( |29| ) then yields the condition 
-R(W 2 -!)(! + b)(W - Rb(l + b)) 



RW(1 + b)+ AW 2 - 2 



(29) 



(30) 



Inserting b and eliminating the unphysical root W = 
—R/2, eq. (|3(]) can be expanded identically into eq. j2^). 

In the mass method of GGB, the condition for circular 
orbits is obtained by equating the ADM and Komar mass 
(as obtained from maximal slicing) 



Madm = M K . 



(31) 



Inserting (|Tl]) and ( |l6| ) yields, after some manipulation 
and elimination of the unphysical root W = 0, again 
the condition (|26|) . Thus we have established that both 
criteria yield the correct condition for circular orbits in 
our model problem. 

Since ( |3l| ) only holds for stationary spacetimes, this 
criterion is closely related to a relativistic virial theorem. 
This relation is also evident from the expansions of the 
ADM and Komar masses to first order in e ~ 1/R ~ v 2 , 



M A 



DM 



1 

2R 



1 



V = 1+U+T 



and 



1 



R 



-if = 1 + 2U + 3T, 



(32) 



(33) 



where U and T are the Newtonian potential and kinetic 
ener gies of the spherical shell. The two expansions fl3^ ) 
and (|33| ) are equal only if the Newtonian virial theorem 
T = -U/2 holds. 

For completeness, we evaluate the relativistic virial 
theorem in spherical symmetry as derived by fijl] 



AnS - 



1 //din 2 l/(ilnV' 2 ^2 
V? [} dr ' ~ 2' dr 



Y)) ^ 6 r 2 dr = 0. 

Computing the above integral in terms of the Komar and 
ADM masses yields 



(W 2 - 1) 

w 



2M 2 



Madm - 2M K + 2R 



%^ = 0, (35) 



which can again be brought into the form (Gf 
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We now briefly discuss the physical implications of the 
condition (Eq). Solving for R we find 



R 



4W 5 - 4W S + W 
6W 4 - 10W 2 + 4 ' 



(36) 



To find a minimum value for the radius of our shell, we 
extremize the above equation with respect to W, which 
yields 

(2W 2 — l)(6W 6 — 21W 4 + IbW 2 — 2) = 0. (37) 

The only physical root (i.e. W real and W > 1) is W = 
1.607, corresponding to -R m i n = 1.532. Expressing this in 
terms of A/adm and circumferential radius Rq we find 



V A/adm / , 



2.506 (equilibrium). 



(38) 



This value should be compared with the Buchdahl limit 
(i? c /Af A DM)min = 9/4 = 2.25 H for static fluid balls 
and (i?c /-Wadm )min = 3 Q for test particles in circular 
orbit in Schwarzschild spacetimes. 

Requiring the particles' orbits to be stable leads to a 
more stringent limit on the compaction, which we find by 
requiring the second derivative of A/adm with respect to 
R to vanish in addition to (27). This yields an equation 
for W with the physical root W — 1.108 corresponding 
to R m i n = 3.053, or 



( Rc 



\M A 



DM 



= 4.265 (stability), 



(39) 



which should be compared with the innermost stable 
circular orbit (Z?c/A/ADM)mm = 6 of test particles in 
Schwarzschild spacetimes. 

To summarize, we construct an analytic solution to 
Einstein's field equations describing a thin spherical shell 
consisting of collisionlcss particles in circular orbits. We 
apply the turning-point criterion ( p7j ) used by CB and 
the mass criterion ( |3~i| ) used by GGB and find that both 
conditions correctly identify circular orbits. The later 
criterion is intimately related to adopting maximal slic- 
ing, which is a natural choice for constructing quasi- 
equilibrium spacetimes (compare flq| ). Our calculation 
illustrates these two criteria in the context of a very 
transparent, analytical framework and provides further 
evidence that the differences between the findings of CB 
and GGB result from the different initial value decom- 
positions. 
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